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NOMENCLATURE 


a particle  radius 

C single  particle  cross  section 

E energy 

f fraction  of  the  snowpack  surface  occupied  by  ice  particles 

H irradiance 

K absorption  coefficient 

N number  of  particles/unit  volume 

n refractive  index 

Q factor  used  to  convert  the  transmission  coefficient  to  the  transmissivity  and  assure  the  con- 
servation of  radiant  energy  at  an  interface 
R reflectivity 

T transmissivity 

X extinction  coefficient 

y radiation  phase  change  produced  by  transmission  across  a boundary 
e emissivity 

r)  radiation  phase  change  produced  by  reflection  at  a boundary 
8 angle  of  incidence  or  coaltitude 
X wavelength 

p density 

a scattering  coefficient 

f skin  depth 

<t>  angle  of  refraction 

& azimuth  angle 

Subscripts 

c cavity,  denotes  the  emissivity  associated  with  the  pores  on  the  snowpack  surface 
i ice,  incident  and  imaginary  used  as  pi(  0;  and  n-x,  respectively 
r real 

s snow,  denotes  properties  associated  with  the  bulk  snowpack 

ab  absorption 

ext  extinction 

sc  scattering 

tr  transmitted 


SNOWPACK  OPTICAL  PROPERTIES 
IN  THE  INFRARED 


Roger  H.  Berger 


INTRODUCTION 


The  interaction  of  electromagnetic  radiation  with  a 
natural  snow  cove-  is  of  considerable  practical  interest 
and  scientific  importance.  The  remote  sensing  com- 
munity requires  an  improved  understanding  of  the  op- 
tical properties  of  snow  in  order  to  correctly  interpret 
the  imagery  of  snow-covered  regions  and  to  distinguish 
between  snowpacks  and  cloud  covers.  Radiant  energy 
transfer  is  important  in  the  thermodynamics  and  meta- 
morphosis of  the  snowpack  and  is  therefore  of  inter- 
est to  those  who  study  the  mechanical  and  hydrologi- 
cal properties  of  the  snowpack.  Therefore,  the  estab- 
lishment of  relationships  between  the  optical  proper- 
ties of  snow  and  frequently  measured  snowpack 
quantities,  such  as  grain  size  and  density,  would  be  of 
considerable  utility.  A model  is  proposed  in  this  re- 
port for  estimating  the  infrared  absorption,  reflection 


and  transmission  coefficients  and  the  emissivity  of  a 
snowpack  from  its  mean  grain  size  and  density. 

As  shown  in  Figure  1,  the  infrared  region  of  the 
electromagnetic  spectrum  extends  from  about  0.8  to 
1000  pm,  bu*  absorption  by  various  atmospheric  con- 
stituents limits  most  infrared  remote  sensing  to  two 
atmospheric  windows  at  3-5  and  8-14  pm.  The  thecy 
of  the  infrared  optical  properties  of  snow  developed  in 
this  report  is  valid  for  wavelengths  between  2 and  20 
pm,  the  region  most  useful  for  remote  sensing. 


THE  MODEL 


Radiation  incident  upon  the  snowpack  surface  under- 
goes either  single  or  multiple  interactions  as  illustrated 
in  Figure  2.  These  multiple  interactions  have  been 
treated  by  Dunkle  and  Gier  (1955)  and  Bergen  (1970) 


Visible 


10  10  10 


|Xfn 


icP  I o'  10° 


1 

' P 

' ‘uHF1  VHF 

X-Roys 

uv  i 
p 

Rodio 

a.  The  electromagnetic  spectrum. 


Expanded 

Area 


b.  Designations  within  the  infrared 
region  of  the  electromagnetic 
spectrum. 


B a?  50t 


Hi 

l/i 

-i20 

l!/> 

U 

iy  I 

rv  ii 

— i ; — i — i i p ■ — i — > i ' — i — i i — i i 

16/im 


c.  Atmospheric  transmission  in  the 
infrared  showing  the  principal  atmos- 
pheric constituents  causing  absorp- 
tion. The  horizontal  arrows  indi- 
cate the  atmospheric  ‘windows" 
through  which  most  remote  sensing 
in  the  infrared  is  done. 


Figure  7.  The  electromagnetic  spectrum. 
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Figure  2.  Interaction  between  electromagnetic  radia- 
tion and,  the  surface  layer  of  ice  particles  in  a snow- 
pack.  The  solid  lines  represent  the  incident  radiation, 
the  dotted  tines  represent  radiation  that  is  specularly 
reflected  either  singly  or  multiply,  and  the  bold  dashed 
lines  represent  refracted  and  absorbed  radiation. 


and  Giddings  and  LaChapelle  (1961 ) as  a radiative 
transfer  or  diffusion  problem  in  a continuum.  If  the 
work  of  these  authors  is  extended  from  the  visible  to 
infrared  wavelengths,  problems  arise  dile  to  the  large 
absorption  coefficient  of  ice  in  the  infrared. 

A much  simpler  approach  to  this  scattering  prob- 
lem and  one  which  is  easily  adaptable  to  the  longer  in- 
frared wavelengths  is  to  look  at  the  problem  in  terms 
of  geometrical  optics  (Bohren  and  Barkstrom  1974). 

A rigorous  treatment  of  the  snowpack  model  proposed 
in  this  report  would  require  Mie  scattering  theory,  but 
the  arbitrary  shapes  of  the  actual  ice  grains  within  the 
snowpack  make  this  refinement  unnecessary. 


SNOWPACK  DESCRIPTION 


The  snowpack  model  proposed  is  a random  distri- 
bution of  spherical  ice  grains  immersed  in  a nonabsorb- 
ing medium.  Below  the  plane  which  defines  the  sur- 
face of  the  snowpack,  the  grain  distribution  is  homo- 
geneous and  isotropic  and  may  be  characterized  by 
the  mean  grain  radius  a,  the  density  of  ice  p-x  and  the 
snow  density  p%.  Since  the  density  of  ice  is  a slowly 
varying  function  of  the  temperature,  p-x  can  be  con- 
sidered a constant.  For  geometrical  optics  to  be 


applicable  to  this  model  two  conditions  have  to  be 
satisficd.  First,  the  particles  must  be  large  compared 
to  the  wavelength  of  the  radiation.  Second,  the  inter- 
granular spacing  must  be  large  compared  to  the  average 
grain  dimensions.  The  first  r-  odition  is  easily  satisfied, 
since  for  most  snows,  the  mean  grain  diameter  is  more 
than  an  order  of  magnitude  larger  than  20  pm,  the 
maximum  wavelength  considered  here.  The  second 
implies  that  the  geometrical  optics  approach  is  valid 
only  for  low-density  snow.  The  range  of  snow  densities 
for  which  this  assumption  is  valid  will  be  examined 
more  closely  in  relation  to  the  results  of  the  theory. 


RADIANT  INTERACTIONS 


If  we  consider  a single  grain  within  the  snowpack,  it 
emits  radiation  isotropically  and  is  immersed  in  the  flux 
of  radiation  from  the  surrounding  grains.  Let  us  con- 
sider this  external  radiation  flux  as  a plane  wave  and 
examine  the  interaction  between  it  and  an  ice  particle. 
This  incident  radiation  of  intensity  H will  either  be 
absorbed  or  scattered  and  these  processes  can  be  des- 
cribed in  terms  of  the  single  particle  cross  sections  c* 
for  absorption  and  Cx  for  scattering.  These  cross  sec- 
tions are  defined  by  the  ratios  of  the  energy  absorbed 
or  scattered  to  the  incident  flux: 


Cab  * Eib/H  for  absorption 


(la) 


and 


C. c 2 £SC/W  f°r  scattering. 


(1b) 


By  the  law  of  conservation  of  energy  these  cross  sections 
are  related  to  the  extinction  cross  section: 


* C,l+Cc 


(2) 


With  the  initial  assumption  of  a random  distribution  of 
grains  spaced  relatively  far  apart,  these  cross  sections 
may  be  related  to  the  bulk  material  absorption,  scatter- 
ing and  extinction  coefficients  by 


*s  = Cext/V- 


(3a) 

(3b) 

(3c) 


The  assumptions  of  large  grain  separation  and  random 
distribution  are  necessary  so  that  the  contributions 
from  the  individual  grains  may  be  added  without  con- 
sidering phase  coherence.  The  number  of  grains  per 
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h'iqure  3.  Ray  paths  illustrating  how  the  Incident  radi- 
ation H interacts  with  an  ice  particle.  The  single  particle 
absorption  cross  section  is  calculated  by  summing  the 
energies  \HT(1-e~Ki),  HTR e~KS(l-e-Kt),...\  which 
are  absorbed  along  successive  ray  paths  [ 1,2  ,...  ] 
of  length  2 within  the  sphere.  The  scattering  cross  sec- 
tion is  calculated  by  summing  the  initially  reflected 
energy  HR  and  the  transmitted  energies  [ HT2e~KS, 
HTlR(<T*®J2,...]. 


AO  ro  8.0  10.0  12.0  14.0  4.0  6.0  8.0  10.0  12.0  14.0 

Wavtlength  (Mm)  Wov.lenglh  f/xm) 

c.  Real  part.  b-  Imaginary  part. 

Figure  4.  Wavelength  dependence  of  the  complex  refractive  index  of  bulk  ice.  (Adapted  from  the  measurements  of 
Schaafand  Williams  1973). 


volume  N may  be  expressed  as  a function  of  the  snow 
density  and  grain  size: 

N = 3ps/4rra3p;.  (4) 

Therefore  the  absorption,  scattering  -nd  extinction 
coefficients  of  the  snow  can  be  expressed  as  functions 
of  the  snow  density  and  grain  size  if  the  single  particle 
cross  sections  for  these  processes  can  be  calculated. 


CALCULATION  OF  CROSS  SECTIONS 

For  calculation  of  the  cross  sections,  consider  a 
plane  wave  incident  upon  a spherical  ice  grain  as  shown 


in  Figure  3 and  on  the  cover.  The  sphere  is  immersed 
in  air  with  a refractive  index  of  1,  and  the  ice  has  a 
complex  refractive  index,  ft=nt-inx.  The  dependence 
of  the  refractive  index  upon  wavelength  is  illustrated 
in  Figure  4.  In  the  succeed’  •’  development  this  wave- 
length dependence  is  implic..  so  that  all  quantities  de- 
pending upon  the  refractive  index  should  be  considered 
spectral  quantities.  The  fractions  of  the  incident  energy 
reflected  and  transmitted  at  each  air/ice  interface  are 
given  by  the  Fresnel  coefficients  R and  T respectively. 
The  irradiance  at  a distance  s along  the  first  ray  path 
inside  the  sphere  (path  1,  Figure  3)  is 

(5) 


»yw*>  ,’vv 


Therefore  the  total  energy  adsorbed  by  the  sphere 
along  path  1 of  length  £ is 

//)ab  = Wr(1-e*K8). 

Similarly  for  path  2 (see  Fig.  3) 

tf2ab  = (Wre-K8)/?(!-e-Ke) 

and  the  total  for  all  internal  transits  is 


(6) 


(7) 


Wab  =Wr(1-e-K8)^(/?c-K(!!i  (8) 


y=i 


7=o 


where  the  sum  is  the  addition  of  the  effects  of  each 
transit.  For  the  total  energy  absorbed  by  the  sphere 
from  the  incident  plane  wave,  'his  expression  is  inte- 
grated over  al!  angles  of  incidence  0 ,. 


/•ir/2  f2n 

^ab=  d0A  < 

J 0 Jo 


d\l/  sinflj  cos0. 


a2HT{\  -e-K8)^"'  [tfe*188]' 


(9) 


where  0 is  the  azimuthal  angle. 

The  total  radiation  scattered  by  the  sphere  may  be 
expressed  similarly: 


/>jr/2  /»2  r 

-=1  dd] 

Jo  Jo 

<72//[/^/?e-K8'V(/?e“K8)i] 


d\li  sinfl:  cosflj 


(10) 


where  the  first  term  is  due  to  the  reflection  of  the  in- 
cident beam  and  the  other  terms  are  due  to  the  radia- 
tion that  is  transmitted  after  1, 2, 3 ...  transits  of  the 
sphere. 

The  reflectivity  and  transmissivity  are  derived  from 
the  Fresnel  formulae  (Born  and  Wolf  1965): 

A Aw 


rA  = R'I*  e,f!i  = (cos0-rtco$0}/{cosG+ 

+ncos^)  (1 1 a) 

tL  = 7f2  e'71  = 2cos0/(cos0+/jcos<j>)  (11b) 


rll  = R'n  e,,!l1  = &COS0~COSJ$)I&'COSQ+ 

1 J J 

i 

it 

jj£ 

+cos$) 

(11c) 

h 

i 

f|(  = Tff  e'7li  = 2cos0/(r/cos  0 +cos$) 

(lid) 

'it 

r 

*} 

■h 

where  tfie  symbols  i and  ||  refer  respectively  to  the  ori- 

>  i 

entation  of  the  electric  vector  perpendicular  and  parallel 
to  the  plane  of  incidence.  In  these  formulae  both  n and 

0 are  complex,  since  they  are  related  oy  Snell’s  law: 

i i 

sin0  = r?sin$! 

(12) 

Substituting  r^cos^u-zV  and  rN/t,.-//?,  into  the  Fresnel 

formulae  and  solving  for  the  reflectivity  and 
sivity  components  yields 

transmis- 

1 

Rl  = [(cos0-i/)2+i'2]/((cos0+t/)2+i/2] 

(13a) 

I 

Tl  = 4co'2O/[(cos0+m)2+v2) 
R..  = 


(13b) 

(13c) 


(nf+n()2  cos20-2\u(nf-nf)+2vnrni\  Cqs0+u2+v2 
(n2+n2)2  cos20+2[u(n2-n2)+2vnrnl\  cosflV+v2 

7*|,  = 4cos20/|(n2+/?2)2  cos20-<  2\u(n2-n2) 

-••2w7r/?i]cos0+tf2+v2}.  (13d) 

In  order  to  satisfy  the  law  of  conservation  of  energy  at 
the  boundary  of  the  ice  particle,  the  fraction  of  the  in- 
cident energy  which  is  reflected  and  the  fraction  trans- 
mitted must  add  up  to  unity.  This  condition,  which 
must  be  satisfied  by  each  polarization  component,  is 
achieved  by  the  addition  of  factors  Q„  and  Ql  : 


Rl+QJl  - 1 and  = 1 


where 


Ql  = uj  cos0 


and 


Ql  = !£/(n2-/7?)+2t'/7rni]/cos«. 


(14) 


(15a) 


(15b) 


For  circularly  polarized  light,  the  polarization  com- 
ponents may  be  combined  to  form  the  total  transmis- 
sivity and  reflectivity: 

T = ZAQt  Tl  +<?., r„)  and  R = V4(^  +*„}-  (1 6) 


& 


/V* 


Combining1/?  cos$  = u-iv  with  Snell’s  law  and  solving  * 
for  u and  v yield: 


u = 707  /?2-/?2 -sin20 


+ ,/(/?2-/??-sin20)2+4/72/?2 
v = O.7O7|-(/?2-/?2-sin20) 

+ V(/?2-/?2-sin20)2+4//2/72 


'h 


J4 


(17a) 


(17b) 


It  will  be  noted  that  neither  eq  13  nor  17  contain 
any  dependence  upon  the  integration  variable  \j/,  so 
that  eq  9 and  10  may  be  integrated.  These  integrals 
may  then  be  simplified  further  by  examining  the  ex- 
ponential terms.  The  absorption  coefficient  of  ice  in 
the  3-  to  1 00 -pm  wavelength  range  is  between  1 00 
and  14,000  cm'1  (Irvine  and  Pollock  1968,  Schaaf 
and  Williams  1973),  and  eq  9 and  10  become 

/-ir/2 

£ab  = 2vo2H  I T(0,nr,nf)s\n0cos0dB 
Jo  (18) 


and 


= 2 iw2H  f 
•n) 


nl  2 


R (0,  nr,  /7j)sin0cos0c/0. 

(19) 


Equation  1 8 gives  the  total  energy  absorbed  by  a 
single  ice  particle  of  radius  a,  immersed  in  a radiation 
field  of  incident  radiant  intensity  H.  The  remaining 
energy  in  the  incident  field  is  scattered  by  reflection 
with  the  energy  distribution  strongly  peaked  in  the 
forward  direction  as  shown  in  Figure  5. 

The  single  particle  cross  sections  for  absorption  and 
scattering  can  now  be  found  by  combining  eq  1 , 1 8 
and  19. 


RESULTS 

As  a result  of  the  large  absorption  coefficient,  all 
the  scattering  of  the  radiation  incident  upon  the  par- 
ticle is  due  to  reflection.  Combining  eq  la,  3a, .4  and 
18  yields  the  absorption  coefficient  for  the  ensemble 
of  ice  particles: 


K.  = 1 .5ps/opj 


r*l  2 


T{0,  nr,  /7j)sin0cos0 dd. 

(20) 


90* 


Figure  5.  The  reflected  energy  distribution 
as  a function  of  the  angle  of  incidence.  The 
energy  has  been  averaged  over  the  wave- 
length interval  for  each  of  the  atmospheric 
windows. 


This  absorption  coefficient  is  dependent  upon  the  ratio 
of  the  snow  density  to  the  particle  size  and  wavelength. 
This  wavelength  dependence  is  illustrated  in  Figure  6 
where  the  quantity  K^/ps  is  plotted  as  a function  of 
wavelength.  Within  the  3-5-  and  8-14-pm  wavelength 
bands  for  typical  snowpack  values  for  the  grain  size 
and  density,  the  absorption  coefficient  is  in  the  range 
of  2.6  to  1 2.7  cm'1 . Under  the  same  circumstances 
the  scattering  due  to  reflection  has  a calculated  range 
of  0.2  to  0.8  cm"1 . 

In  addition  to  the  absorption  and  reflection  due  to 
the  refractive  index  of  the  particles,  the  diffraction  ef- 
fects must  be  included  when  considering  forward 
scattering.  All  the  energy  which  is  incident  upon  the 
particle  is  removed  from  the  incident  wave  either  by 
absorption  or  reflection,  giving  an  effective  cross  sec- 
tion equal  to  the  projected  area  of  the  sphere,  to2. 

By  Babir.ct’s  principle  the  cross  section  for  light  dif- 
fracted by  a large  particle,  i.e.  1<<  x < 2to/X,  is  equal 
to  its  projected  area.  Therefore  the  total  extinction 
cross  section  is  twice  the  geometrical  cross  section: 


^ext  = ■ 


(21) 


WSlHVt:- 
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Figure  6.  The  normalized  absorp  tion  coefficient  as  a function 
of  wavelength.  The  arrows  indicate  the  two  atmospheric  win- 
dows, 3-5  and  8-14  pm. 


The  extinction  coefficient  for  snow  is 


EMISSIVITY  CALCULATION 


*s  = /Cs+ots  (22 ) 

where  the  subsctipt  ts  of  the  scattering  coefficient  a 
denotes  total  scattering  effects  of  the  particle.  Follow- 
ing the  notation  of  Bohren  and  Barkstrom  (1974)  an 
effective  extinction  coefficient  may  be  defined  by 

X*  - Ks+or+ot  (23) 

where  the  a and  t stand  for  the  reflected  and  trans- 
mitted parts  of  the  scattering  coefficients.  Since  in 
the  derivation  of  eq  19  from  eq  10  it  was  seen  that 
the  transmitted  radiation  was  negligible,  ot  is  essen- 
tially zero.  From  Babinet’s  principle  and  eq  3c  the 
effective  extinction  coefficient  is 


X*  = O.lSpJap-,  (24) 

and  by  writing  the  scattering  and  absorption  coeffici- 
ent in  terms  of  this  parameter  we  have 

or  = 0.01 3X*  and  Ks  = 0.99 X*  (25) 

which  shows  that  Q9%  of  the  effective  extinction  is 
due  to  absorption. 

Defining  an  extinction  and  an  absorption  depth  as 
the  reciprocal  of  the  coefficient,  these  depths  are  re- 
spectively 

rext  = 50  mm  and  tKs  = 3.8  mm. 

These  values  are  calculated  for  a typical  aged  snow- 
pack  with  density  of  350  kg/m3,  mean  crystal  size  of 
2 mm  and  ice  density  of  91 7 kg/m3. 
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This  absorption  length  is  of  the  order  of  a few  par- 
ticle diameters  and  indicates  that  no  significant  amount 
of  radiation  at  these  wavelengths  penetrates  beyond 
the  first  few  ice  particles.  The  snovvpack  emissivity 
can  therefore  be  found  by  considering  only  the  radia- 
tive interaction  with  these  surface  layers.  To  account 
for  all  the  radiation  that  is  absorbed  within  the  surface 
layers,  three  absorption  processes  have  to  be  considered, 
These  processes  are  1 ) absorption  at  initial  incidence 
with  the  ice  particles,  2)  absorption  after  one  or  more 
reflections  within  the  snowpack,  and  3)  cavity  aosorp- 
tion  in  the  pore  spaces  within  the  snowpack.  These 
processes  are  expressed  by  the  first  two  terms,  the 
third  and  fourth  terms,  and  the  last  term,  respectively, 
in  the  following  equation  which  is  derived  using  eq  1 8 
and  19: 


XT XlZ 

T(d,n',n-J  sin  0 cos  6 dB 

/Sin'1  (r/o) 

r(0,r?r,/7j)sin  0 cos  6 dd 


f*l  2 

+2na2HfA  I R(0,nr,n;)sin0cos0d0 

Jetp) 


„* -i 


/•sin'  M 

+2ftr2HPBi  R[0,nvn^s\n6cosB  dB 

Jo\P) 


+ec[l-7rAr2/3(o2+r2)] 


(26) 


Figure  7.  The  spectral  emissivity  of  snow  illustrating  the  de- 
pendence  upc " the  density,  ps  in  g/cm3. 


where 


/•if/2 

lH\  [ 

*'o 


1/4=2  TO2///  \R(9wfrn*jttjnx)\ 


sin0cos0d0 


and 


]/B  = 2irr2H 


i 


sin '{r fa) 


[R{0,nt,t.t)+T(0,  nvn$] 


sin0cos0</0 


and  where  f is  the  fraction  of  the  surface  area  occu- 
pied by  the  projection  of  the  surface  layer  of  particles, 
and  P is  the  fraction  of  the  surface  area  occupied  by 
the  portion  of  the  projection  of  the  subsurface  par- 
ticles which  is  not  obscured  by  the  surface  layer  of 
particles.  The  two  functions  d(p)  and  0‘(p)  define  the 
minimum  angle  of  incidence  at  which  total  absorption 
takes  place  within  .the  snowpack  structure.  ec  is  the 
apparent  emissivity  of  a pore  which  opens  onto  the 
surface  of  the  snowpack. 

In  order  to  evaluate  these  integrals  it  is  necessary 
to  make  a further  assumption  about  the  snowpack. 
This  assumption  is  that  the  average  properties  of  the 
snowpack  may  be  approximated  by  a regular  array  of 
particles.  Since  most  of  the  radiation  incident  upon 
an  individual  ice  particle  is  absorbed  and  the  absorp- 
tion length  is  short,  the  effects  of  phase  coherence 
may  be  neglected.  Under  these  conditions  0(p)  is  a 
nearly  linear  function  of  density  and  varies  between 
50.33°  and  51 .36°  over  the  density  range  of  1 70  to 
481  kg/m3.  The  second  function  0(p)  is  a constant 


22.5°  over  this  range.  The  two  projected  areas  as  func- 
tions of  the  density  become 

f=  (Mf>//3  (27) 

f>  = 7r[1/2-v/2(/Wps),/3+(Afps)2/3]  (28) 

where  M = 3/4ttP|.  Substituting  eq  27  and  28  into  eq 
26,  the  expression  for  the  emissivity  becomes 

e = Tr(Afps)2/3{Min20(p)+(1  l^2-(Mps)m  j 2 1 

(Mps)2/3-sin2p/8| 

+27r(Afps)24f  n T(0,nr,n)sin0cos9  dO 


fe(p) 

1 


] 


+ i 7"(0,/7r,/j[)sin0cos0</0 


h{l-Jt(1/2-V2  (Afps)1/3+2(Mps)2'3]|< 


(29) 


The  pore  emissivity  ec  is  0.95,  based  on  the  calcula- 
tions of  Cbandos  and  Chandos  (1974).  The  evaluation 
of  eq  29  for  the  3-14-pm  spectral  region  is  shown  in 
Figure  7.  It  can  be  seen  from  this  figure  that  the  de- 
pendence of  emissivity  upon  the  snow  density  is  very 
weak.  Likewise,  the  emissivity  is  not  a function  of 
the  mean  grain  size,  a. 

There  have  been  very  few  emissivity  measurements 
on  snow  in  the  2-14-pm  wavelength  range  and  there 
are  no  published  results  known  to  this  author  which 
report  its  spectral  emission  or  reflectance  and  characterize 
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the  snow  by  grain  size  and  density.  Values  of  the 
emissivity  appearing  in  the  literature  range  between 
0.81  to  0.99,  with  the  higher  value  being  most  com- 
mon, but  generally  no  information  is  included  on  -he 
spectral  bandwidth  or  the  snow  properties.  The  most 
complete  documentation  is  furnished  by  Dunkle  and 
Gier  (1955)  who  measured  the  emissivity  of  two 
natural  and  two  artificial  snows  and  obtained  values 
of  0.82,  0.89,  0.81  and  0.95.  The  lower  values  apply 
to  the  finer-grained  snow  of  each  variety.  The  differ- 
ence in  emissivity  between  the  fine-  and  coarse-grained 
snow  was  attributed  to  either  particle  size  or  surface 
roughness,  but  the  change  in  density  might  have  also 
been  the  cause  since  no  density  measurements  were 
made.  These  results  are  not  sufficiently  detailed  for  a 
definitive  comparison  with  the  theory  developed  in 
this  report. 

For  a definitive  test  of  this  theory  it  is  necessary  to 
have  the  snow  samples  well  documented  so  that  the 
effect  of  grain  size  and  density  variations  ma/  be 
known.  The  wavelength  band  of  the  measurement 
must  also  be  limited  to  the  region  of  validity  of  the 
theory.  This  is  important  because  inclusion  of  radia- 
tion of  wavelengths  shorter  than  2.5  /xm  will  reduce 
the  measured  emissivity  since  reflectance  measure- 
ments in  this  region  indicate  that  the  emissivity  de- 
creases rapidly.  A series  of  emissivity  measurements 
which  can  be  compared  with  this  theory  is  in  progress. 
Preliminary  results  indicate,  in  agreement  with  the 
theory,  that  the  emissivity  does  not  depend  upon  the 
grain  size,  but  the  dependence  on  density  is  not  dearly 
defined.  These  results  will  be  presented  in  a subsequent 
report. 


DISCUSSION  OF  RESULTS 


The  theory  used  to  derive  the  optical  constants  of  a 
snowpack  contains  no  adjustable  parameters  and  is 
based  only  upon  the  refractive  index  of  pure  bulk  ice 
and  the  density  and  average  grain  size  of  the  snow- 
pack.  In  order  to  evaluate  the  validity  of  the  results 
derived  above,  the  implications  of  the  initial  assump- 
tion of  spherical  grains  and  large  interparticle  spacing 
must  be  examined. 

If  an  arbitrarily  shaped  particle  is  compared  with  a 
spherical  particle  of  the  same  volume,  the  sphere  has  a 
smaller  surface  area.  Chylek  (1977)  has  shown  that 
the  average  geometrical  cross  section  of  the  nonspheri- 
cal  particle  is  always  larger  than  that  of  the  correspond- 
ing spherical  particle.  It  is  uncertain  what  the  net  ef- 
fect of  this  on  the  emissivity  would  be,  since  tha  last 
term  in  eq  29  would  be  decreased  while  the  other 
terms  would  be  increased.  Since  the  absorption  coef- 


ficient of  ice  is  large,  the  absorption  length  within  an 
ice  particle  is  much  smaller  than  the  particle  dimen- 
sions and  therefore  the  absorption  coefficient  is  inde- 
pendent of  the  particle  shape. 

The  limitations  of  the  theory  due  to  the  assumption 
of  large  interparticle  spacing  are  more  difficult  to  assess. 
The  development  of  the  bulk  material  constants  from 
the  particle  interaction  cross  sections  is  dependent  upon 
the  particle  spacing  being  large  enough  so  that  phase  co- 
herence effects  may  be  neglected.  Under  these  condi- 
tions the  contribution  of  each  particle  may  be  added  to 
get  the  scattered  field  of  the  ensemble.  The  random 
distribution  of  the  particle  sizes  (Bohren  and  Barkstrom 
1974,  Hodkinson  and  Greenleaves  1962)  and  the  ran- 
dom shape  and  orientation  (Bohren  and  Barkstrom 
1974)  of  the  particles  within  the  snowpack  may  average 
out  the  phase  coherence  effects  as  the  density  of  the 
snow  increases.  Experiments  with  pigments  (Blevin 
and  Brown  1961)  show  that  the  reflectance  of  an  en- 
semble is  independent  of  particle  concentration  ever 
a wide  range  and  then  usually  decreases  at  high  con- 
centrations, indicating  that  the  theory  developed  for 
snow  may  be  valid  for  densities  up  to  400  kg/m3. 

Therefore  the  absorption  and  scattering  coefficients 
and  the  emissivity  values  calculated  using  the  theory 
developed  here  should  be  useful  for  densities  between 
170  and  400  kg/m3.  The  theory  also  shows  that  the 
absorption  and  scattering  coefficients  are  linearly  re- 
lated to  the  snowpack  density  and  inversely  related  to 
the  grain  size.  The  emissivity  is  independent  of  grain 
size  and  exhibits  only  a weak  dependence  upon  the 
snowpack  density. 
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A facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 
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